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Abstract 

Weak interactions within a nucleus generate a nuclear spin dependent, parity violating electro- 
magnetic moment, the anapole moment. We analyze a method to measure the nuclear anapole 
moment through the electric dipole transition it induces between hyperfine states of the ground 
level. The method requires tight confinement of the atoms to position them at the anti-node of 
a standing wave Fabry Perot cavity driving the anapole-induced micro-wave El transition. We 
explore the necessary limits in the number of atoms, excitation fields, trap type, interrogation 
method, and systematic tests necessary for such measurements in francium, the heaviest alkali. 

PACS numbers: 32.70.Cs, 31.15.Ar, 32.80.Pj 
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I. INTRODUCTION 



Zel'dovich postulated in 1957 that the weak interactions between nucleons would gen- 
erate a parity violating, time reversal conserving moment called the anapole moment 
Flambaum and Khriplovich calculated the effect it would have in atoms 2|]. Experiments 
in thallium gave a limit for its value {s], and it was measured for the first time with an 
accuracy of 14% through the hyperfine dependence of atomic parity nonconservation (PNC) 



m cesium 



We present in this paper a measurement strategy of the nuclear anapole moment by di- 
rect excitation of the microwave electric dipole (El) transition between the ground hyperfine 
levels in a chain of isotopes of an alkali atom. Alkali atoms are the best understood quan- 
titatively in their electronic properties associates with PNC. The precision of the Cs PNC 
experiments has required more detailed studies of the nuclear structure Q|. Measurements 
over a chain of isotopes offer the advantage that they can focus on the differences appearing 
as the number of neutrons changes. This task has been accomplished well by theory (see for 
example Ref. (t]) for the hyperfine anomaly measurements in Fr. 

Current plans at the Isotope Separator and Accelerator (ISAC) at TRIUMF, in Vancouver 
Canada, should provide access to all the neutron deficient long lived isotopes of Fr with 
lifetimes above 30 s and to a similar number of the neutron rich isotopes, a sufficient variety 
to give a difference in number of neutrons of more than 10. The expected production rates 
should be at least two orders of magnitude larger than those obtained at Stony Brook, the 
leading place for study of Fr [Sj. 

A measurement of the anapole moment in a chain of isotopes will provide information 
about neutral weak currents in the nucleus. The measurements can also give information 
o. the nudea. structure a.d its changes a. the „u„,be. of .eut.o,. ...eases flQ. 

The El transition between hyperfine levels is parity forbidden, but becomes allowed by 



the anapole induced mixing o 



suggested in the past [U 



levels of opposite 



13, Q; 



15 



16|, 
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mrity. The general approach has been 



18l |. We would place many atoms inside 



a microwave Fabry-Perot cavity and hold them in a blue detuned dipole trap. The atoms 
would interact with the microwave field and with a Raman field generated by a pair of 
laser beams, in the presence of a static magnetic field. We would confine the atoms to the 
node (anti-node) of the magnetic (electric) microwave field to drive only an El transition 
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between hyperfine levels. The atoms would start in the lower hyperfine level, with the signal 
proportional to the population of atoms in the upper hyperfine level after the excitation. 



The interference with a Raman transition would give a signal linear in the E 



transition. 



Recent work related to time-reversal invariance tests in atomic traps |19l . |20[ |. points to 
the many potential advantages of combining traps with tests of fundamental symmetries, 
but also highlights the potential systematic errors present in such measurements, making a 
careful evaluation of the method prior to its implementation necessary. We focus our study 



p..a., on isotope, of .ano.,n, t.e heaviest of the a.ka>, ato„. fl, . an optica, dipo.e 

trap, where the effect is expected to be large. 

The organization of the paper is as follows: Section II gives the theoretical background for 
the nuclear anapole moment, section III explains the proposed measurement method, section 
IV presents an analysis of noise sources and systematic effects, and section V contains the 
conclusions. 



II. THEORETICAL BACKGROUND 

The exchange of weak neutral currents between electrons and nucleons constitute the 
main source of parity violating atomic transitions. The currents are of two kinds, depending 
on whether the electron or the nucleon enters as the axial vector current. The corresponding 
terms in the Hamiltonian differ on their dependence on the nuclear spin. The part indepen- 
dent of the nuclear spin is generally the dominant contribution in atomic PNC. This is not 
the case for the present work, where we consider transitions between hyperfine levels of the 



ground state, and the contribution from the nuclear spin independent part is zero 2l| . The 



Hamiltonian for the spin dependent part in the shell model with a single valence nucleon of 



unpaired spin is given by 



22| 



where G = 10~^ m~^ is the Fermi constant, m^ is the proton mass, K = (/+l/2)(— 1)'^+^/^^', 
/ is the nucleon orbital angular momentum, I is the nuclear spin, a are Dirac matrices, and Ki 
is the interaction constant, with i = p,n for a proton or a neutron. The terms proportional 
to the anomalous magnetic moment of the nucleons and the electrons have been neglected. 



The interaction constant is given by 22| 



_ K - 1/2 / + 1 

with K2,p = —i^2,n = —1.25(1 — 4sin^ ^^H^)/2, corresponding to the tree level approximation, 
with sin^ 9\Y ~ 0.23 the Weinberg angle. Equation [2] has two corrections, the effective 
constant of the anapole moment, and kq^ that is generated by the nuclear spin indepen- 
dent part of the electron nucleon interaction together with the hyperfine interaction. The 
three parts of this interaction constant can be traced to different ways in which the weak 
interacting vector boson Z° appears in the Feynman diagrams. The first one (the anapole) 
correspond to vertex corrections due to weak hadronic interactions on the nuclear side of 
the electromagnetic interaction coupled to the electron through a virtual photon. The sec- 
ond one takes the direct effect of a exchange between the electron vector current and 
the nuclear axial current. The last one is the simultaneous exchange of a and a photon 



modifying the hyperfine interaction. Flambaum and Murray showed that 22 1 



9 «^2/3 
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where a is the fine structure constant, yUj and /ijv are the magnetic moments of the external 
nucleon and of the nucleus respectively, ro = 1.2 fm, A is the atomic mass number, Qy^ is 
the weak charge, and gi gives the strength of the weak nucleon-nucleus potential with (^p ~ 4 
for a proton and 0.2 < < \ for a neutron 23|]. The anapole moment is the dominant 
contribution to the interaction in heavy atoms, for example, i^a,p/i^Qw — 1^ for 209pj,_ 
will assume from now on that k,- = k„ . 



A. The anapole moment 

The anapole moment of a nucleus is a parity non-conserving (PNC), time reversal con- 
serving moment that arises from weak interactions between the nucleons (see the review by 



Haxton and Wieman 2J].) It can be detected in a PNC electron- nucleus interaction, and 
reveals itself in the spin dependent part of the PNC interaction. Wood et al. {4, 5] measured 
the anapole moment of ^^^Cs by extracting the dependence of atomic PNC on the hyperfine 
levels involved. 
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The anapole moment classically is defined by (see for example 25l |) 



TT / d\r^3{v), (4) 



with J the electromagnetic current density 26|. The nuclear anapole moment in francium 



arises mainly from the weak interaction between the valence nucleons and the core. Flam- 
baum, Khriplovich and Sushkov estimate the anapole moment from Eq. H] for a single 
valence nucleon to be 

G K 

where j is the nucleon angular momentum and e the electron charge. The calculation 
assumes a homogeneous nuclear density, and a core with zero angular momentum, leaving 
the valence nucleon carrying all the angular momentum. 

The measurement of the anapole moment gives information on the weak nucleon-nucleon 
interactions. A measurement of the anapole moment in a chain of isotopes would provide a 
separation of the anapole moment due to the valence proton or neutron. 

B. Calculations of the anapole moment of francium isotopes 

We use Eqs. [3] and [5] to estimate the anapole moment of five light francium isotopes with 
radioactive lifetimes longer than one minute 7|]. The unpaired valence proton generates the 
anapole moment in even-neutron isotopes, whereas in the odd-neutron isotopes both the 
unpaired valence proton and neutron participate. Francium has an unpaired hg/2 proton for 
all the isotopes and a /5/2 neutron for the odd-neutron isotopes around ^^°Fr. The protonic 
and neutronic contributions add vectorially to generate the anapole moment: 

C-j; ■ I + j„ G (/ + 1/2) 

with Cf"'ji the anapole moment for a single valence nucleon i (proton or neutron) as given 
by Eq. [5](jp = 9/2, j„ = 5/2). Equation [6] defines the coupling strength of the total anapole 
moment (/€«) resulting from adding the valence proton_and neutron. Figure [1] shows the 
predicted values of K,a for a string of francium isotopes |lS] using = 
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FIG. 1: Anapole moment effective constant for different isotopes of francium. 
C. Perturbation theory 

The anapole moment induces a small mixing of electronic states of opposite parity. The 
effect of the anapole moment Hamiltonian on the ground state hyperfine levels according to 
first order perturbation theory is 



\7F^) = \sFm) + Y: ^-^^^^4^^4^\pF'm'), 



F'm' 



Ep — Eg 



(7) 



where Ep, and Eg are the energies of the p, and s states respectively, and 



Ha = |e|cK-a5(r), (8) 

is the anapole moment Hamiltonian from Eq. [H with a the anapole moment from Eq. O 
The matrix element in Eq. [7] gives 23 1 

. ^Z'^R 27 + 1 (/ + l/2)KaRy 



{pF'm'l Ha \sFm) 



{g^Qpy/^ 3 + 
x{F{F + 1) - /(/ + 1) - 3/4)5^,^,5„,„,, 



(9) 



with ^ = Gw?^o?l\pl'n = 3.651 x 10"^'', me the electron mass, Z the atomic num- 
ber, Qs and Qp the effective principal quantum number for the s and p electronic states, 
7 = + 1/2)^ — Z^a^, J the electron total angular momentum, and Ry the Rydberg. 
The relativistic enhancement factor R is given by 



i? = 4(ao/2Zro)'-'Vr2(27 + l), 



(10) 



with ao the Bohr radius, and Tq = TqA^^^. 



The anapole moment mixes only states with the same F and m, and the mixing grows 
as TF'/'^R. For the ^ospj^ ground state, we obtain 



\sFm) = \sFm) - i 5.9 x lO'^^Ka 

x(F(F + 1) -25.5)|pFm). (11) 

The mixing coefficient is imaginary due to time reversal symmetry. In practice, the mixing 
would be measured through the El transition amplitude Aei (Eq. [16]) it induces between 



two hyperfine levels. The effect in francium is 11 times larger than in cesium 27 1. 



III. PROPOSED MEASUREMENT STRATEGY 



High efficiency magneto-optical traps (MOT) for francium atoms on line with an accel- 



erator have been demonstrated 28j. Their performance and reliability matches the needs 
of the current proposed measurement strategy. Atoms captured on a ffist trap would then 
be transferred to a second MOT in a separated chamber. We would load the atoms into a 
dipole trap located at the electric field anti-node of a standing wave in a microwave Fabry- 
Perot cavity. We would optically pump them into a single Zeeman sublevel, and prepare 
a coherent superposition of the hyperfine ground levels with a Raman pulse of amplitude 
Afj and duration t/j. Simultaneously we would drive the El transition of amplitude Aei 
with the cavity microwave field, and measure the population in the upper hyperfine level 
(normalized to the total number of atoms [A^]) using a cycling transition. The population 
in the upper hyperfine level at the end of each sequence would be 

" ATI |2 AT ■ 2 f i^R^ ^El)tR\ . s 

^± = N\ce\ =Nsm ( — 1, (12) 

where Ce is the upper hyperfine level amplitude. The sign depends on the handedness of 
the coordinate system defined by the external fields, as explained in the next section. The 
signal for the measurement: 



S = E_ 



. (ARtn\ . (AEitR 
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N.i.[^)(^], (13) 
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FIG. 2: Schematic setup of the proposed apparatus. The microwave cavity axis is along the y-axis. 
The microwave electric field inside the cavity oscillates along the x-axis. The two Raman laser 
beams are polarized along the x-axis and z-axis, respectively. The microwave magnetic field and 
the static magnetic field are both directed along the z-axis. A dipole trap (not shown) holds the 
atoms at the origin that coincides with an anti-node of the microwave electric field. 

would be the difference between populations in the upper hyperfine level for both handed- 
ness. The last step assumes a small Aei, the quantity proportional to the anapole moment 
constant Kq. 

A. Apparatus setup 

Figure [2] shows a diagram of the proposed apparatus. The atoms would be placed inside 
a microwave Fabry-Perot cavity at the electric field anti-node, confined in a blue detuned 
dipole trap to a volume with 10 fim length along the cavity axis, and a 1 mm diameter in 
the radial dimension. Observation of the electric dipole (El) microwave transitions would be 
done through an interference method and extraction of the signal would require repeating 
the excitation varying the coordinate system. 

Preparation of the atoms in a particular Zeeman sublevel of the lower hyperfine level 
|Fi,mi) in an applied static magnetic field B = Bqz would be necessary. A resonant 
standing-wave microwave electric field Ei{t) = E cos{2TTh'mt + 4') cos(/cml/)x would excite the 
atoms to a particular Zeeman sublevel in the upper hyperfine level \F2, m2). The microwave 
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magnetic field M would be aligned along B, and it is 7r/2 out of phase (for a perfect standing 
wave) with E so that M(t) = Msin(27rz/mt + ip) sm{kmy)z, with M = E in cgs units. 

The Raman transition would include two plane- wave optical fields, E/ji(t) = 
Ejii cos(ti;i(;t + 0_r)x and E/j2(i) = -E/?2 cos{{ljji + 2TxVm)t + phase locked to the mi- 

crowave field. The Raman carrier frequency uo^ would be tuned sufficiently far from optical 
resonance that only the vector part of the Raman transition amplitude (V oc iEiri x '£1112) 



would be non-negligible 



29| : that is, we ignore the tensor part of the Raman amplitude. 



B. Observable and reversals 



The various electric and magnetic fields of the apparatus would define a coordinate system 
related to the measured rate 'E±. The transition rate S-t depends on three vectors: The 
polarization of the El transition (E), the polarization of the Raman transition (V), and the 
static magnetic field B that provides an axis for the spins of the nuclei. We combine these 
three vectors to produce the time reversal preserving pseudo scalar 2(E x (E^^xE/j2) ■ B), 
proportional to the measured quantity. 

A single reversal of any of the fields in the above pseudo scalar changes the sign of the 
interference term of Sj-. We then would have the following reversals: 

1. - Magnetic field reversal {(3 reversal). 

2. - A shift of vr in the relative phase between the El and the Raman fields (s reversal). 
The Zeeman sublevels reverse with the magnetic field. The state preparation has to be 

inverted in order to reach the correct Zeeman sublevel, meaning that light goes into 
cr^ and vice versa. The magnitude of the static magnetic field and the microwave cavity 
frequency remain unchanged for this reversal. 



C. Apparatus requirements 

1. Magnetic field 

We would drive El transitions between two particular Zeeman sublevels, |Fi,mi) — ^ 
1^2,7712) in different hyperfine levels of the ground state. While the frequencies of the excit- 
ing fields can be well controlled, the energy difference of the Zeeman states is determined 
primarily by the static magnetic field. 
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TABLE I: Parameters of the five relevant francium isotopes: Spin, hyperfine splitting (Hfs) of the 



7si/2 state |3d . 



31[ , Zeeman sublevels mi , m2 , and their energy separation at the static magnetic 



field Bq used in the proposed measurement. 



Isotope Spin Hfs(MHz) 


mi m2 


i?o (Gauss) 


Z^m(Mhz) 


208 


7 


49880.3 


0.5 1.5 


2386.5 


49433 


209 


9/2 


43033.5 


-1 


1553.0 


42816 


210 


6 


46768.2 


0.5 1.5 


2586.4 


46208 


211 


9/2 


43569.5 


-1 


1572.3 


43349 


212 


5 


49853.1 


0.5 1.5 


3265.7 


49015 



The experimental design should minimize the sensitivity to magnetic field fluctuations. 
The energy difference between two levels passes through a minimum at the static magnetic 
field Bq, and depends quadratically on the magnetic field around that point. We would use 
the Zeeman sublevels that give the smallest quadratic dependence. Table |T] lists the Zeeman 
sublevels and magnetic fields selected for different francium isotopes. The experiment would 
work between the |Fi,mi) and \F2,'m2) levels and also between the |Fi,m2) and |F2,mi) 
levels, interchanging mi and m2. The operating point of the static magnetic field and the 
frequency of the microwave cavity would have to be corrected slightly because of the nuclear 
spin contribution. The state preparation would also change to start in the appropriate level. 
The change of mi (m2) for m2 (mi) does not work as a reversal because of the difference in 
transition amplitude, but it can still be useful as a consistency check. 

The frequency for the F = 4, m = to the F = 5, m = — 1 transition in 209pj,^ expanded 
around the critical field Bq = 1553 Gauss, is 



Um = 42.816 X 10^ + 90(5 - BofRz, (14) 

with B in Gauss. Control of the magnetic field to 0.06 Gauss (three parts in 10^) reduces 
the frequency noise due to magnetic field fluctuations down to ~ 0.3 Hz. 

The experiment would take place in a large magnetic field whereas the state preparation 
and detection occur in a small magnetic field. The transition between both regimes should be 
done adiabatically. The time scale is determined by the precession time in a small magnetic 
field, resulting in a magnetic field ramp duration of hundreds of microseconds. 

10 



2. The microwave cavity 



The francium hyperfine separation requires a Fabry-Perot microwave cavity operating at 
around 45 GHz (wavelength Am ~ 0.66 cm) in a Fabry-Perot configuration; for example a 
cavity with a mirror separation of d ~ 20Xm ~ 13 cm and a mirror radius of = 3.5 cm. 
These parameters combine to minimize diffraction losses as the Fresnel number > 1, 
where = r'^/Xmd |32[ |. 

The quality factor (Q) of the cavity is 

Q = |, (15) 

where <^ is the skin depth and is equal to ^J^jujii^a with /xq the magnetic constant and a 
the conductivity (5.8 x 10^ fi^^m^^ for copper at room temperature). The conductivity 
limited quality factor is Q = 1.9 x 10^. It is possible to couple 58 mW into the cavity with 
current available technology, which would give an electric field of 476 V/cm to drive the El 
transition. 

The El transition amplitude for 209pj^ between the initial hyperfine level (i) F = 4, m = 
to the final hyperfine level (/) F = 5, m = — 1 with a static magnetic field of 1553 Gauss 
(see Table H]) is 



(/|-eE-rK)/n 
O.Oli 



E 



476 V/cm _ 



0.45 



rad/s. 



(16) 



A more accurate result can be obtained with the use of many-body perturbation theory 
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271, 



34|. 



A 1 cm cavity waist would cover the atoms in the 1 mm diameter 10 /im length trap, 
and radius of curvature of Rm = 9.9 cm for the cavity mirrors ensure a stable cavity, since 
(1 — {d/2Rm)Y < 1. The curvature of the wave fronts could create a gradient of polarization 
of the microwave field smaller than 3 x 10^^ rad cm~^ over the volume of the trap. We show 
later that this rotation is within acceptable ranges. 

The field inside the cavity can be decomposed into a standing wave and a travelling 
wave. The presence of the travelling wave generates Ml transitions despite the location 
of the atoms at the node of the standing wave magnetic field. Significant reduction of the 
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amplitude of intra-cavity travelling waves comes with a symmetrical arrangement of identical 
antennas, one on each minor. Antennas give a higl, coupling efficiency into the cav,ty 

as compared to a slit or a grating ^]. The electric field inside the cavity is given by 



X [Eiti (e*'"'""' - r2e'^"'''-e~'^"''') 

+E2t2 (e"*^'-" - ne^^'-'^e^'^'"")] , (17) 

where r is the reflectivity, t the transmissivity, k is the wave-vector of the microwave field, 
d the separation between the mirrors, and the sub indices 1 and 2 refer to the two mirrors. 
The first (second) term is the field generated by antenna 1 (2). The expression is the sum of 
two waves, one travelling to the right and the other to the left. The difference in amplitude 
between these two contributions results in a travelling wave. The ratio of travelling to 
standing wave assuming a symmetrical cavity, that is ri = r2 = r and ti = ^2 = t, is 

nr/s =(j + ^^^) (^(1 + KAd) , (18) 
with d the phase mismatch from both antennas and Aci the deviation of the cavity mirrors 
separation from the ideal position. Assuming d = Q and control of the amplitude from 
each antenna to 1%, the position of the mirrors to 0.1 /im and taking 1 — r = 3.6 x 10~^ 
(consistent with the Q factor computed above), we obtain 1Zt/s = (3 + 9i) x 10"^. 



3. Dipole trap 

We choose a far-detuned dipole trap to contain the atoms for the duration of the mea- 
surement since the perturbations introduced by it are small and measurable. A variety of 
different geometries have been proposed over the years. These include red-detuned traps 



371, 



38| for 



based on focused beams, and blue-detuned traps with hollow beams (see Refs. 
reviews of recent work.) 

The trap would confine the atoms within 10 /im around the microwave electric field 
anti-node and 1 mm diameter in the radial dimension. The region of confinement would 
be smaller than the microwave wavelength (Lamb-Dicke regime), so Doppler broadening 
becomes negligible. 
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The AC Stark shift (AE"), which produces the restoring force of the dipole trap, displaces 
the two hyperfine levels of ground state in the same direction but not by the same amount. 
The differential shift changes the resonant frequency for the cavity-driven El transition 
used in the anapole moment measurement. The change in the hyperfine separation for a 
detuning {6 = w — We) larger than the hyperfine splitting (Ahfs) is approximately equal 
to (A//i?s/(5)Ai? [39|. The shift reduces considerably using a blue detuned far off resonance 
trap (FORT) at 532 nm. 

The dipole trap in combination with the cavity field may generate a multi- 
photon transition. There are four vectors available for that transition: EId, 
MId the dipole trap electric and magnetic fields, E the microwave electric field 
and B the static magnetic field. The parity and time reversal conserving ob- 
servables created with combinations of the above vectors that produce a res- 
onant transition ((EId • E)(M1d ■ B), (EId ■ B)(M1d ■ E), (EId x E) ■ (MId x B), 
(EId X B) ■ (MId x E), (EId x MId) ■ (E x B), and 2(E1d x E) ■ MId), give a negli- 
gible contribution if the trap laser propagates along B. 



4. Ml transition 

The dominant transition between the two hyperfine states is a magnetic dipole Ml tran- 
sition. The magnetic component of the microwave field could drive Ml transitions. A 
microwave magnetic field polarized along the x axis would have the same signature as a 
parity violating signal. The Ml transition amplitude {Ami) between the levels of interest is 
given by 



{f\{-e/2me){J + S)-M\t)/h 
M 



7.8 X 10' 



1.6 Gauss 



rad/s. 



(19) 



for the maximum expected microwave magnetic field in the Fabry-Perot cavity. The ratio 
of the El transition (Eq. [T6|) to the Ml transition is ~ 1 x 10^^. The success 

of the measurement depends on reducing and understanding this transition. We propose to 
suppress it in three ways. 

First (see Fig. [3[a]), we would place the atoms at the magnetic field node (electric field 
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anti-node) of the microwave cavity. The magnitude of the microwave magnetic field at the 
edges of the atomic trap is reduced by a factor K = sin(27r(ij/Am), with dt = 10 fim the 
length of the trap along the cavity axis. The reduction factor at 45 GHz is K = 4.8 x 10"^. 

Second (see Fig. E^b]), we would direct the polarization of the Ml field to be along the 
2-axis (Fig. [2]). The non-resonant Ml transitions in this case would be of the type Am = 0. 
The static magnetic field (Bq) would split the Zeeman sublevels of the two hyperfine levels, 
and the microwave field would be resonant for the \Am\ = 1 El transitions (the microwave 
electric field would be polarized along the x axis.) The alignment imperfections give a 
suppression factor equal to sin(0) ~ ~ 10"^ rad, the angle of the microwave magnetic 
field polarization with respect to the z axis. 




, dipole trap 

Brf / potential 



FIG. 3: Suppression mechanisms of the Ml transition, (a) Trapped atoms would sit at the magnetic 
field node, where the magnetic field is zero, (b) Schematic of the ground hyperfine levels showing 
a \Am\ = 1 transition such as the one for the anapole moment and the Am = out of resonance 
transition such as that induced by the Ml field. The level spacing as well as the spin do not 
correspond to any particular atom, (c) Trapped atoms would oscillate around the microwave 
magnetic field node and would sample a zero time-averaged magnetic field. 

Third (see Fig. EJ^c)), the atoms in the dipole trap would oscillate around the microwave 
magnetic field node. An atom crossing the node would see a microwave magnetic field 



14 



pointing in the opposite direction. The change in position effectively would flip the phase of 
the magnetic field that the atom sees, and would reverse the evolution generated by the Ml 
transition. The dynamical suppression only takes place if the frequency of oscillation {() 
of the atoms inside the trap is larger than the Rabi frequency of the Ml transition and is 
given by {1/^/N)^Imi/C- The frequency of oscillation along the cavity axis for the proposed 
geometry would be C/Stt ~ 300 Hz. 

Taken together, the three suppression mechanisms would reduce the expected Ml transi- 
tion amplitude to Amis/^ = 1-9 x 10^^ rad/s for 10^ atoms. This is 500 times smaller than 
the amplitude for the El transition. 



D. Signal to noise ratio 

The magnitude of the signal from Eq. [13] reaches a maximum for a Raman transition 
amplitude of = {2n + l)vr/2 with = 1. The measurement of the upper hyperfine 
state population collapses the state of each atom into one of the two hyperfine levels. The 
collapse distributes the atoms binomially between the two hyperfine levels and leads to an 



uncertainty in the measured excited state fraction called projection noise Afp 
projection noise is given by 



3- 



The 



Afp = v/iV|Ce|2(l-|Ce|2). (20) 

The projection noise vanishes when all the atoms are in one of the hyperfine levels, but in 
those cases the noise is dominated by other sources, such as the photon shot noise. 
The signal to noise ratio for a projection noise limited measurement is 

4- = 2^^Vn. (21) 

Up h ^ ' 

Taking Aei from Eq. [161 ^i? = 1 s, and integrating over 10^ cycles, we would reach a 3% 
measurement with only 300 atoms. 

The high-efficiency MOT that we developed at Stony Brook, with production rates around 



10^ s~^, captures in excess of 10^ francium atoms 28|]. We expect to trap 10^ atoms after 
transferring them to an ultra-high vacuum environment. In this case, Eq. [21] predicts a 
signal to noise ratio of 20 in 1 s. Higher francium production rates could be obtained at 
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other facilities, such as ISAC at TRIUMF, where an actinite target could deliver in excess 
of 10^ atoms per second of a single isotope. 

While measurements in francium benefit from a large A^i, large atomic samples of other 
alkalis are easily prepared. We could obtain the same signal to noise ratio in a cesium sample 
with 100 times more atoms and the same strength-driving field. While the fundamental 
signal to noise ratio indicates the inherent trade-offs between different alkali species, technical 
noise, specific to the instruments dedicated to the measurement, must also be considered. 
For a discussion of technical noise in the cesium PNC Boulder experiment see Ref. j4|. 

IV. NOISE AND SYSTEMATIC EFFECTS 

The measurement of the anapole moment would come from determining the population 
transferred from the lower to the upper hyperfine level by the application of the Raman and 
microwave fields. Both of these fields (or any other stray field) are characterized by a field 
amplitude, frequency (or detuning), and interaction time. The total transition amplitude 
for a common detuning (S) and interaction time {tn) is: 

A = {Am + Aeu + A,) + i{AR2 + Aeu + ^2), (22) 

where Aiji j?2 are the real and imaginary components of the Raman amplitude, Aeii,ei2 the 
corresponding for the El transition amplitude and Ai^2 are the real and imaginary parts of 
any other transition present such as an Ml transition. 

Table [TTl shows the phase of the transitions for given field polarizations, with their trans- 
formation under magnetic field reversal assuming all the excitation fields are in phase. We 
control the phase difference between the Raman field and the cavity El field. Varying 
ip introduces an additional factor of e^"^ on the El transition amplitude while the Raman 
transition remains unchanged. The standing wave Ml field inside of the cavity is 90° out of 
phase with the El field, which gives a factor of ie^"^ for the Ml transition. If instead the Ml 
field corresponds to a traveling wave, then it is in phase with the El field. 

The Raman field would be polarized along the y axis so that Am = A^y, and the El 
transition polarized along the x axis so that Aeu = ^Aeix (or = vr/2). These two 
amplitudes would interfere since both are in phase and only one (the El) changes sign under 
magnetic field reversal as shown in Table [TTl Expanding Eq. [12] for large A^y compared to 
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TABLE II: Phase {P = A/\A\) of the relevant transition amplitudes for the initial state Fi = 4, 
nil and final state F2 = 5, m2 and polarized along the specified axis. For this table all the fields 
have the same phase (equal to 0). P/ja; represents the Raman transition with one vector along the 
y axis and the other along the z axis, such that their cross product points along the x axis. /3 
represents the static magnetic field reversal together with a sign change on the Zeeman sublevel 
m. 

Reversal ml m2 Peix Pmix Pmiv Prx Prv 
Normal 0-1 i 1 i i 1 
/3 1 -i -1 i -i 1 

the detuning and other amphtudes, we obtain 
with, 

Aef={iAEx + Ai + —— + -—{AR^ + A2f]. (24) 

V -^^Ry ^Ry J 

Aef contains the signal (Aex) and noise {Ai, A2, A^ and 5) terms. We can use this expres- 
sion to set limits in the different experimental parameters and identify the corresponding 
observable. Expanding the last term in Eq. [23] for small Ir gives 

^^i^y^e/. (25) 

The first term in A^f is proportional to iAjiyAEix, which corresponds to the PNC signal 
^(Ex (E^iXEk2)-B). 

The amplitudes of interest are the Raman amplitudes Aiix,Ry, the El amplitude Aeix, a 
Ml transition that is in phase with the El field AMix,Miy and an Ml transition that is 7r/2 
out of phase with the El field AMox,Moy As an example, if the standing wave magnetic field 
inside of the cavity is tilted towards the x axis it generates an amplitude Amox since this 
field is out of phase with the El field. The Ml amplitudes are included in Eq. |2l] as Ai or 
A2 depending on their phase relation to Ajiy. 

The relevant values for the relative phase {ip) between the El and the Raman transition 
are multiples of it/2. First we study the case with ip = 0,7t. This does not correspond to 
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the PNC measurement since the El and Raman transitions are out of phase and do not 
interfere. The signal obtained with this configuration is still helpful in the evaluation of 
unwanted contributions. We can rewrite Af,f from Eq. [21] using Table [TTl and ignoring the 
detuning (5) as 



+ s{iAfiyAMoy — '^iAfixAMox) + (^{ — '^iAMiyAMox) 

+sf5{ARyAMix - 2AM^yAn^)] , (26) 

with s = 1, — 1 when = 0,7i respectively and (3 = 1, —1 depending if we have the normal 
experiment or we apply a magnetic field reversal. With = 7r/2, —tt/2 instead we get 



^e/ — . [(^Mia; ^tloy) 
^Ry 

+ s{iAjiyAMiy + ^iAji^AMix) — Pi—^iAMoy^Mix) 

+sp{2AR^AMoy - AnyAuox + lARyAEix)] , (27) 

where now s = 1,-1 when ip = 7r/2, —tt/2 respectively. This corresponds to the experimen- 
tal condition for the PNC measurement. The PNC signal is contained in the last term, and 
it changes sign under both reversals. Equations [26] and [27] show how reversals can be used 
to isolate the PNC signal. 

We divide the analysis of the different experimental parameters into three parts: Sys- 
tematic effects that include terms that mimic the PNC signal and that are contained in the 
last parenthesis of Eq. [271 line broadening mechanisms, which contain all other terms and 
that average to zero after an infinite number of cycles, and calibration errors that modify 
the value of the extracted constants on the PNC signal. 



A. Line broadening mechanisms 

We start with terms that do not change under both reversals. They include the detuning 
term from Eq. [2l]and all the terms in Eq. [27] except for the last parenthesis. We present the 
requirements to achieve a precision of 3% in the measurement after 10^ repetitions. Each 
noise amplitude has to be controlled to 3Aei 
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We could reduce the effect of some noise terms by increasing A^y (see Eq. |2l]). We would 
take Apiy to be exactly equal to a {2n + l)7r/2 pulse, and include any deviation from this 
value into Ai. We would control the Raman pulse to 0.025% in one second with shot noise 
limited detection. This would limit the maximum value for the Raman pulse to Apty/h = 121 
rad/s or n = 38. We now proceed to analyze the spurious terms in Eq. [27]that contaminate 
the signal. 



1. h'6'/2ARy 

The detuning can have its origin in poor frequency control on the microwave or Raman 
fields, or changes in the external fields that shift the energy levels. The detuning would have 
to be controlled to 6 = 2.7 rad/s. The required accuracy for the microwave field frequency 
is one part in 10^^. 

Control of the static magnetic field Bq to a fractional stability of 5 x 10~^ would keep 
the detuning under control. 

The presence of an Ml transition produces an AC shift of the levels. The value of the 
maximum shift is ~3 mHz, which is negligible. 

The atoms in the trap occupy different vibrational levels. Transitions between different 
vibrational levels are suppressed for a sufficiently far detuned trap. Each vibrational level 
has slightly different resonance frequency that leads to broadening of the signal and loss of 
coherence. 

Coherence times as long as 4.4 s have been measured for atoms in a blue detuned trap [41^. 
The main source for decoherence was the distribution of Stark shifts felt by the atoms. We 



expect a coherence time 16 times smaller in francium than in Ref. j41| using a laser at 532 
nm because of the difference in hyperfine splitting and detuning. The dephasing grows slowly 
in time and can be reversed with the use of an "echo" technique. The atoms would spend 
approximately half of the time in each hyperfine level with a Raman transition amplitude 
Ar = {2n + l)vr/2 for large n. It is necessary in that case to keep the coherence for a time 
approximately equal to t/j/n, with tji the duration of the experiment. We would need a 
coherence time of 26 ms for n = 38 to have an interaction time of 1 s. This is below the 
expected 300 ms coherence time. 

The average differential Stark shift seen by the atoms would be approximately equal 
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to kT{AHFs/S)/h = 6.3 Hz. The effect of the time varying detuning generated by the 
oscillations in the trap is similar to a steady state detuning of the same magnitude, and can 
be compensated by adjusting the microwave frequency. We must control the power of the 
trap laser to 7%. 

2. A\JAny 

This term appears due to a bad polarization alignment of the Raman field. Control of 
the polarization of the Raman field to one part in 10^ would be necessary to suppress this 
term. 

These terms are multiplied by a small number and their contribution becomes negligible. 
For example, Ajij.AMix/ ^Ry has the small factor Aji^^/Ajiy appearing due to polarization 
misalignment in the Raman beams. 

4- AnyAMiy/ARy = AMiy 

This is the dominant term that depends on the Ml transition. The Ml field appears due 
to imperfections in the microwave cavity field that create a traveling wave component that 
may be in or out of phase with the El transition. 

Eq. [18] gives the amplitude of the traveling wave expected in our setup. The traveling 
wave is polarized along the z axis, so we can include the polarization suppression factor of 
10^'^. Combining these two numbers with the amplitude for the Ml transition we get an 
amplitude of 0.25y4Ei out of phase with the El transition, and an in phase amplitude of 
0.75Aiji. 

The relative phase between both antennas [d) can be adjusted by minimizing the Ml con- 
tribution when the static magnetic field (B) is tilted slightly. The antennas phase mismatch 
contribution remains controlled for d < 0.01 rad. 
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B. Systematic effects 

The systematic effects include the terms in the last parenthesis in Eq. [271 They change 
sign under both s and (3 reversals just as the PNC signal. The constraints for these terms 
are stronger since they do not average to zero. Their contribution must be below 0.03^^1 
to reach a 3% measurement. We proceed to analyze each one of these terms. 

1- ARxAMoy/ARy 

This term appears because of a combination of misalignment of the Raman field and 
misalignment of the microwave field or imperfections in the microwave cavity. It corresponds 
to the observable M x (E^-^xEr2) ■ B. This term is reduced by the Raman misalignment 
i^Rx/^Ry) and its contribution would become negligible. 

2- AnyAMox/ARy = AmOX 

This term has the same origin as the previous one, but its contribution is considerably 
larger since it is not suppressed by the Raman misalignment. It gives the limiting factor 
in the precision of the measurement and its control depends completely on the suppression 
mechanisms. 

The cavity mirrors may have some birefringence, which generate a microwave magnetic 
field X-axis component. The microwaves make roughly 1000 reflections in the cavity. We 
need a polarization rotation smaller than 10~^ rad or a rotation per reflection smaller than 
10^^ rad to keep the Ml suppression unchanged. The constraint for a 3% measurement is 
14 times smaller. 

The atomic sample would have to be precisely held at the node of the microwave magnetic 
field. The maximum displacement we can tolerate is 3 x lO^^^m for a 3% measurement. 

C. Calibration errors and requirements on theoretical calculations 

The PNC signal (Eq. [13]) would give directly the Aei amplitude since the uncertainty 
in the Raman amplitude is negligible. A^i is the product of the microwave electric field 
and the matrix element. The microwave electric field amplitude has to be known to 3%. 
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The electric field could be measured by tilting the magnetic field and inducing an Ml 



transition. The extraction of information a 
measurement requires theoretical input 24 



30ut the weak interaction from an experimental 



42l | . The quality of the electronic wave functions 



is the most important. The accuracy of the matrix elements has to be comparable to 
that of the experiment. The effective constant of the anapole moment is obtained after 
subtracting the other two contributions to Hi (Eq. [2]). Johnson et al. show that the other 
contributions for the case of Fr amount to a few percent [27!]. The anapole moment of the 
even-neutron isotopes comes only from the unpaired proton, while the odd-neutron isotopes 
contain contributions from the unpaired proton and neutron. A measurement of the anapole 
moment to better than 10% would give an initial separation of both contributions 22|. 



D. Other sources of fluctuations 



The microwave magnetic field would generate transitions to other levels of the type 
Am = 0, which are non-resonant at the proposed magnetic field (detuning ~ 0.4 GHz). 
Nevertheless, these transitions will have to be taken into account in a detailed analysis of 
the data. 

Stray electric fields produce Stark induced transitions that mimic the PNC signal. A 
stray electric field of 13 V/cm in the z direction would generate a transition amplitude 
equal to the parity violating signal. Stray fields large enough to be a problem are unlikely 
to occur and can be ignored [29 1. 

Gradients induce higher order multipole transitions, such as an E2 transition. Fortu- 
nately, these higher order transitions between the two hyperfine ground levels are strongly 
suppressed. Table IIIII summarizes the results of the analysis of noise and systematic effects. 



V. CONCLUSION 



The anapole moment provides a unique probe of weak hadronic interactions. In particu- 
lar it is sensitive to weak long-range meson exchange interactions, and consequently allows 
a measurement of weak neutral currents in the nucleus. This is not the case in high-energy 
experiments where the weak contribution must be separated from the strong and electro- 
magnetic contributions that are much larger. We have presented the analysis of a proposed 
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TABLE III: Fractional stability required for a 3% measurement. The observable associated with 

each constraint is also included. 

Observable Constraint Set value Stability 



ARyAEl 


Microwave amplitude 


476 V/cm 


0.03 


^Ry^Ry 


Raman amplitude 


121 rad/s 


2.5 X 10-^ 




Microwave frequency 


45 GHz 


10-11 




Dipole trap Stark shift 


6.3 Hz 


0.07 




DC Magnetic field 


1500 Gauss 


4.7 X 10-5 




Raman polarization 


rad 


10-3 rad 


ARyAMiy 


Mirror separation 


13 cm 


7.7 X 10-^ 




Antenna power 


57 mW 


0.02 




Antenna phase 


rad 


0.01 rad 


ARyAMoX 


Mirror birefringence 


rad 1 


X 10-^ rad 




Trap displacement 


m 3 X 10-11 m 



measurement strategy of the nuclear spin dependent part of the PNC interaction, dominated 
by the anapole moment. While the proposed measurement method can be extended to other 
alkali atoms, a series of measurements in a chain of francium isotopes allows the separation 
of the proton and neutron contributions to the anapole moment. 

As noted by Fortson et al. j^, [lo| studies of atomic parity non conservation give in- 
formation on the nuclear physics. The nuclear weak interaction at low energies is often 
parameterized by a series of coupling constants, either with a meson exchange formalism, 



he so called DDH parametrization 



431] , or more recently with effective field theories (EFT) 
44| . A program of measurements of the anapole moment in a chain of francium isotopes will 
contribute significantly to constraint some of the DDH parameters, which together with the 
EFT program will provide a model independent input for theoretical analysis of low energy 
weak interaction constants. It is important to note that the measurement of the anapole 
moment of an even and an odd isotope of francium give almost orthogonal bands in the me- 
son coupling parameter space. This is subject to the assumption that the anapole moment 
is carried mainly by the last nucleons [22.], but as shown by the measurements of the hy- 
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perfine anomaly |7|], this is a reasonable assumption. These measurements will significantly 
contribute to deepen our understanding of the nuclear structure. 
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